Counter Chemotactic Flow in Quasi-One-Dimensional Path 
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Quasi-one-dimensional bidirectional particle flow including the effect of chemotaxis is investigated 
through a modification of the John-Schadschneider-Chowdhury-Nishinari model. Specifically, we 
permit multiple lanes to be shared by both direct ionally traveling particles. The relation between 
particle density and flux is studied for several evaporation rates of pheromone, and the following 
results are obtained: i) in the low-particle-density range, the flux is enlarged by pheromone if 
the pheromone evaporation rate is sufficiently low, ii) in the high particle-density range, the flux 
is largest at a reasonably high evaporation rate and, iii) if the evaporation rate is at the level 
intermediate between the above two cases, the flux is kept small in the entire range of particle 
densities. The mechanism of these behaviors is investigated by observing the spatial-temporal 
evolution of particles and the average cluster size in the system. 



Recently, models of various particle flows such as traf- 
fic flow, granular flow, and pedestrian flow have been 
investigated Ijj-Q. In general, these systems often experi- 
ence a 'jammed state' with an increase in element density. 
On the other hand, a population of ants taking on chemo- 
taxis with pheromone seems to be able to avoid such state 
even when their density is highjol-fTsj. Thus, such biolog- 
ical aspects seem to provide some hints to eliminate jam 
in particle flows. 

In this study, we investigate the possible effects of the in- 
teraction among particles through pheromone on particle 
flow. We consider a model of bidirectional (i.e., counter) 
particle flow in which the particle deposit pheromone and 
determine their motions according to the pheromone field. 
The model is given as a modification of a cell model of the 
bidirectional ant traffic introduced by John et al.[9|] (here- 
after, the JSCN model). Practically, to realize a wider 
class of flows rather than restricting what to traffic flow, 
we simply extend the JSCN model as follows: I) The sys- 
tem consists of two lanes where the particles traveling in 
both directions share both lanes and are permitted to 
move within a lane and to shift to the other, namely, a 
particle can swerve on encountering the oppositely travel- 
ing another particle. II) The particles momentarily can go 
backward contrary to their traveling direction according 
to the rule described below, where the traveling direction 
is defined as the direction of destination cells at the right 
(or left) edge of the system. 

Details of the model are given as follows. The system 
consists of two lanes and each lane involves L cells. From 
a cell at the left edge to that at the right edge in each 
lane, we give the cell index, 0, 1, L — 1. To realize the 
excluded volume effect, only one particle is permitted to be 
contained per cell. The particle deposits pheromone on its 
occupying cell. The number of pheromone is given zero or 
one which increases with the secretion and decreases with 
evaporation. 

In each sub step, one cell is randomly selected among 
2 X L cells. If no particle is found but a pheromone is 
found in the selected cell, the pheromone in that cell is 
removed with a probability /, whose quantity is defined 
as the evaporation rate. On the other hand, if the se- 
lected cell is occupied by a particle, one candidate cell for 



its next position is determined among neighboring unoc- 
cupied cells. The selection among these neighboring cells 
is given in the following designations order: 1) forward 
cell (traveling directional cell), 2) diagonally forward cell, 
3) side cell, 4) diagonally backward cell, and 5) backward 
cell, as shown in Fig. [H After the candidate cell for the 
next position is determined, it is checked if a pheromone 
is left or not (one or zero) in the cell. If a pheromone is 
left, the particle moves with the probability Q, and in the 
opposite case, with the probability q where Q > q to re- 
flect chemotaxis. After moving in this step, a pheromone 
is put in its newly located cell, and a pheromone in the 
leaving cell is removed with the probability /. When the 
particle reaches the destination cells of the right (or left) 
edge of the system, it changing its traveling direction to 
the opposite jl4|. The same sub steps are repeated 2 x L 
times. All of the above processes constitutes the unit time 
step. 

In the following, we show typical behaviors of this 
model. As the initial condition, TV particles are set ran- 
domly in the field. The parameters are Q = 0.95, q = 0.25, 
and L 500. 
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FIG. 1: (Color online) Ordering for the selection of the moving 
direction of a particle at each sub- step. 



First, we determine the relationship between the par- 
ticle density p = N/(2L)^ the evaporation rate of the 
pheromone /, and the flux J. Here, J is given by the 
long-term average (averaged over time steps ~ 10^) of 
the number of particles reaching the destination cells per 
unit time step. Figure [2] shows the relationship among 
them. As roughly guessed from this figure, several char- 
acteristic types of fundamental diagram are realized in dif- 
ferent ranges of evaporation rates. Here, specific values. 



namely, / = 0, / = 0.01, / = 0.3, and / = 1.0, are chosen 
and the fundamental diagrams realized at each respective 
/ are shown in Fig. [3l 

For / = 0, J linearly increases with p as J = Qp within 
the density region p < 0.07, which reflects the acceleration 
by a pheromone through the relation Q > q. However, J 
falls below the linear relation over a critical density p = 
0.07, which is much smaller than the lowest jam density 
given in the JSCN model because of the lane, sharing rule 
by both directionally traveling particles that makes jam 
more tenacious. 

At higher evaporation rates, / = 0.3 and / = 1, in 
the low-density region, J Hnearly increases with J = qp^ 
which means flux is kept less than that in the above case 
of / = 0. In the case of / = 0.3, however, the linearly 
increasing region extends up to p = 0.09, this is, larger 
flux than the case of / = is realized in a finite range 
beyond p = 0.07 as seen in Fig. 3. In this way, in a 
certain density range beyond p = 0.07, a larger flux is 
realized with a reasonably higher evaporation rate than 
that in the no-evaporation case. 

At an intermediate evaporation rate, / = 0.01, J is 
barely larger than J = in a small-p region. At ap- 
proximately p = 0.055, J begins to decrease, but the shift 
from the increasing region to the decreasing region is not 
so sharp. In this case, J is kept low throughout the entire 
range of densities. 
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FIG. 2: (Color online) Relationship among density (horizontal 
axis), evaporation rate of pheromone (vertical axis) and flux. 
A brighter tone means a higher flux. 



Next, to understand the mechanism giving these char- 
acteristic fundamental diagrams, we determine the spa- 
tiotemporal evolution of particles and the average cluster 
size for each situation. To this purpose, we define the 
term 'cluster', here, as a group of more than one particles, 
which is inseparable by more than one contiguous empty 
columns. Here, a 'column' is a set of cells having identical 
cell indexes. Furthermore, 'entangled cluster' is defined as 
a cluster consisting of oppositely traveling particles that 
stays at almost the same location. 

First, we see the cases of no evaporation, / = 0. At 
a density p = 0.04, particles move without forming large 
clusters, as shown in Figs. |4] and O and each particle 
moves at a velocity close to Q. In Fig. [5l the red/blue 
dots indicate that only right/left directionally traveling 
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FIG. 3: (Color online) J — p relationship; for / = 0(+)5 
O.Ol(O), 0.3(A) and l.O(n). 
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FIG. 4: (Color online) Average cluster size as a function of 
/>, for the evaporation rate: / = 0(+)? 0-01(0)5 0.3(A) and 

i.o(n). 

particle(s) is (are) located in the corresponding columns, 
and green dots indicate that both directional particles are 
contained in those columns. Figures [6]l9] are plotted in 
the same manner. As seen in Fig. 5, when two oppo- 
sitely traveling particles or small clusters encounter, they 
smoothly weave through the other particle or cluster and, 
accordingly, avoid jam formation. At a higher density 
p = 0.1, through the encounter of two oppositely traveling 
clusters larger than those in the above case [Fig. |4], they 
easily form an entangled cluster. The size of this entangle 
cluster grows with time through frequent encounters with 
other incoming clusters [Fig. [6l^a)], and this cluster fi- 
nally absorbs most of the particles in the system, thereby 
blocking flow [Fig. [6fb)]. Then J drastically decreases 
to almost zero [Fig. [3]. Between these two densities, at 
p = 0.075, the system alternately stays in jammed and 
smoothly flowing states [Fig. El^a)]. Then, the average 
flow J takes an intermediate value. 

The above result is unlike the dynamics as reported in 
the JSCN model in which the jammed state containing 
most particles in the system is not dissolved once formed. 
This characteristic dynamics in the present model is real- 
ized because particles are permitted to move against the 
traveling direction according to their temporal situations. 
By combining this backward motion of particle with the 
stochastic selection of a cells performing each sub step, 
the system behaves like a double- well potential system in 
which a particle goes back and forth between two stable 
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FIG. 5: (Color online) Typical spatiotemporal behaviors of 
particles for p = 0.04 and f = 0. The broken line indicates the 
inclination Q meaning the trajectory of free particles in a field 
filled with a pheromone. 
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FIG. 7: (Color online) Spatial-temporal behaviors of parti- 
cles: (a) transitions between two different typical states, (b) 
the initiation of entangled cluster, and (c) the termination of 
entangled cluster, where (b) and (c) are picked up from (a) for 
p = 0.075 and / = 0. 



FIG. 6: (Color online) Typical spatiotemporal behaviors of 
particles in (a) transient with the initiation of entangled clus- 
ter, and (b) the steady state for /> = 0.1 and / = 0. 

states, though they are not exactly the same. 

Below, we discuss the cases of a high evaporation rate 
/ = 0.3. At a density p = 0.06, particles behave simi- 
larly to those in the case of / = and p = 0.04. Namely, 
they move without making large clusters, as seen in Fig. 
m Each particle moves at an average velocity close to q 
[Fig. [8] (a)] which reflects the high evaporation rate of a 
pheromone. When two oppositely traveling particles (or 
small clusters) encounter, they smoothly weave through 
each other to avoid jam formation. As shown in Fig. HI 
at this evaporation rate, the clusters are kept small until 
p = 0.09, accompanied by a smoothly flowing state and a 
linearly increasing part in the fundamental diagram in Fig. 
[3l The reason for the smooth flow being sustained up to a 
higher density than that in the case of / = is yet to be 
studied. At a high density p = 0.12, particles behave sim- 
ilarly to those in the case of / = and p = 0.1. Namely, 
through the encounter of oppositely traveling two clusters, 
they easily form an entangled cluster that grows with time 
through frequent encounters with other incoming of other 
clusters, thereby finally absorbing most of particles in the 
system to block flow. Therefore, J becomes almost zero 
[Fig. [Hb)]. Between these two densities, at p = 0.1, the 
system alternately takes one of the above mentioned two 
states. 

Then, we consider the case of an intermediate evapora- 
tion rate, / = 0.01. At p = 0.04, around which J almost 
linearly increases with p, particles move as medium-sized 




FIG. 8: (Color online) Typical spatiotemporal behaviors of 
particles for (a) p = 0.06 and (b) p = 0.12 with / = 0.3. The 
broken line indicates the inclination q, which is the trajectory 
of free particles in a field without a pheromone. 

clusters (Fig. |4]) with a velocity between q and Q. When 
two oppositely traveling clusters encounter, they tempo- 
rally form an entangled cluster; however, it is immediately 
resolved without growing to include most particles in the 
system. Before long, a new entangled cluster is formed in 
the system to repeat the same process. At a critical den- 
sity p = 0.055, J decreases. However, unlike those in the 
cases of / = 0.0 and / = 0.3, the spatiotemporal behavior 
of the system does not markedly change at the critical den- 
sity, nor the average cluster size discontinuously changes. 
As seen in Fig. O even beyond the critical density, en- 
tangled clusters are not sustained as in the cases of / = 
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FIG. 9: (Color online) Typical spatiotemporal behaviors of 
particles for (a) p — 0.04 and / = 0.01 and (b) p — 0.07 and 
/ = 0.01. 

and / = 0.3; instead, their formation and dissolution are 
repeated in the system, although the average cluster size 
gradually increases with density. Hence, flux is kept low 



but not zero, as seen in Fig. [3l 

In this study, the dynamical and statistical aspects of 
two lanes bi-directional particle flow including the effect 
of chemotaxis were investigated. We found the following 
relationships among particle density, the evaporation rate 
of pheromones, and flux: I) In the low-density region, the 
flux of particles is largest if the evaporation rate is suffi- 
ciently low; II) in the higher-density region, a reasonably 
high evaporation rate enhances flux; and III) flux is kept 
small at all density region if evaporation rate is interme- 
diate between those in cases I and II. The analysis of the 
balance (or the breakdown of the balance) between the 
growth rate and dissolution rate of entangle clusters seems 
to be the key to further investigations, as has already been 
mentioned in a previous study by Kunwar et al.[ll] 
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